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The longitudinal response of single semiflexible polymers to sudden changes in externally applied 
forces is known to be controlled by the propagation and relaxation of backbone tension. Under 
many experimental circumstances, realized, e.g., in nano-fluidic devices or in polymeric networks or 
solutions, these polymers are effectively confined in a channel- or tube-like geometry. By means of 
heuristic scaling laws and rigorous analytical theory, we analyze the tension dynamics of confined 
semiflexible polymers for various generic experimental setups. It turns out that in contrast to the 
well-known linear response, the influence of confinement on the non-linear dynamics can largely be 
described as that of an effective prestress. We also study the free relaxation of an initially confined 
chain, finding a surprising superlinear ~ i 9//8 growth law for the change in end-to-end distance at 
short times. 
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I. INTRODUCTION 

Semiflexible polymers such as actin or microtubules 
play essential roles for the elastic behavior of cellular- 
structures [TJ [2]- Their mechanical properties are de- 
termined by their relatively high bending stiffness, char- 
acterized through a compared to the contour length L 
large persistence length t p . Many of the static and dy- 
namic properties of single filaments have been investi- 
gated in the last decades, and many of the results carry 
over to more complex structures such as polymer net- 
works [TJ [3H5] and solutions [6TTL4] , the latter of which 
are theoretically described through single filaments con- 
fined to an effective tube. In recent years, much experi- 
mental research is also focused on single molecule stud- 
ies, by means of powerful new nano-technological meth- 
ods, where filaments are routinely confined to micro- or 
nanometer-sized channels [T5T425| . 

Theoretical understanding of confined semiflexible 
polymers has mainly been gained in terms of equilib- 
rium statistics [26H34j and some dynamical properties 
in microcapillary flows [35 and on the linear response 
level [5Bl438j . However, for larger forces such as those 
commonly applied in single- molecule experiments, it is 
important to account for the extremely high stretching 
stiffness of most semiflexible polymers, which as an effec- 
tive inextensibility has not only significant consequences 
for the static force-extension relation but also heavily 
influences dynamical properties. The responsible non- 
linearities are especially apparent when externally ap- 
plied forces are suddenly changed, for in this case the 
backbone tension, which prevents stretching of monomer 
bonds, shows a highly nontrivial spatial and temporal 
dependence. Previous approaches to capture this phe- 
nomenon were based on heuristic scaling arguments and 
effective theories [39ti42j . These early results were later 
confirmed and generalized by means of a systematic for- 
malism based on a multiple scale perturbation theory 



[4"3Tl45| . Recent extensions of these results include pre- 
stressed filaments [46], transverse forces [47], extensible 
backbones 48J, and oscillatory forces [3H] • A common 
observation of these studies is that the essential nonlin- 
earities contained in the nontrivial tension profile lead 
to a mixing of short- and long-wavelength fluctuations, 
such that experimentally changing the equilibrium mode 
spectrum of contour fluctuations can have far more dras- 
tic consequences as usually suspected, for instance on the 
scenario-specific relaxation dynamics observed for poly- 
mers that were prepared in initially straight conforma- 
tions by different means [50] . 

It is well known that confinement leads to a changed 
fluctuation spectrum, because contour undulations with 
a wavelength longer than Odijk's deflection length ~ 
£p^ 3 D 2 / 3 (D denoting the channel diameter) [ST] are 
strongly suppressed. Under many circumstances, this 
merely leads to a renormalization of some characteris- 
tic length and time scales [3S1 GUI IS2]- However, in 
the light of the above we expect substantial changes for 
the dynamics of the tension along the contour. Because 
of the widespread presence of confinement in common 
single-molecule experiments, we analyze the effects of 
confinement on the longitudinal response of semiflexi- 
ble polymers in three generic scenarios (cf. Fig. [I]). 
First, we consider the sudden stretching of a semiflexible 
chain in cylindrical confinement upon applying a con- 
stant "pulling" force at its ends. Next, we discuss the 
inverse "release" scenario, where this stretching force is 
suddenly removed. Compared to the corresponding un- 
confined scenarios analyzed previously, we find in both 
cases dynamic signatures strongly reminiscent of a pre- 
stress. In order to clearly elucidate the specific dynamical 
consequences of the different contour statistics of a chain 
confined to a channel, we finally address the idealized 
scenario where this constraint is suddenly removed. We 
will refer to this latter scenario as "free relaxation from 
confinement" (FRC). 
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The paper is organized as follows. In the next sec- 
tion [n] we present the basic mathematical formulation of 
the problem and derive the equations of motion. The 
linearized version of these equations of motion are inves- 



tigated on a scaling level in section III There the physi- 
cal principles behind all subsequently treated relaxation 
processes are illuminated and scaling laws governing the 
longitudinal dynamics are derived. Section |IV| contains 
the analytical core of this work. After a brief survey of 
the ideas that lead Hallatschek et al. [43J [44] to a re- 
fined formulation of the relaxation dynamics of semiflex- 
ible polymers, we adjust this description to the particu- 
lar scenarios outlined above and investigate the tension 
dynamics within the various asymptotic limits. The fol- 
lowing section [V] builds upon these analytical results to 
discuss the longitudinal dynamics of the polymer's pro- 



jected end-to-end distance. Section VI quantifies our an- 
alytical findings in terms of numerical estimates for time, 
force, and length scales relevant in experiments. Finally 
section IVHI summarizes the main results of this work. 
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FIG. 1: Illustration of the three scenarios investigated in this 
work. Top (pulling): Both ends of an initially tension free 
backbone, which is equilibrated in a channel of constant di- 
ameter, are pulled apart by an externally applied force f. Mid- 
dle (release): The confined polymer is equilibrated under an 
externally applied prestretching force f, which is suddenly re- 
leased. Bottom (free relaxation from confinement): The con- 
tour is equilibrated in confinement and relaxes in free space. 
In all cases, the sudden change in external conditions prop- 
agates along the contour within growing boundary layers of 
size £|| (t). 



II. EQUATIONS OF MOTION 

To determine the polymer's equations of motion, we 
resort to the wormlike chain model |53j in its contin- 
uous version [53]. According to this model the poly- 
mer is described by means of an inextensible continu- 
ous space curve r(s) and a bending Hamiltonian "Hbnd = 
k/2 J dsr" 2 introducing an energetic penalty propor- 
tional to the chain's (squared) curvature r" 2 to account 
for bending rigidity (primes indicate derivatives with re- 
spect to the arc length coordinate s). In three dimen- 
sions the modulus k is related to the persistence length 
£ p , which characterizes the decay of tangent-tangent cor- 
relations along the contour, via k — ksT£ p . Confinement 
in a hard-walled channel will be modeled by means of an 
additional harmonic potential H po t = 7/2 j dsr± pun- 
ishing fluctuations rj_ away from the symmetry axis of 
the channel. The relation between the effective confine- 
ment strength 7 and the diameter D of the corresponding 
hard-walled channel is discussed in section VI (cf. Ref. 



|28j). Further introducing the backbone tension f{s,t) 
as a Lagrange multiplier function ensuring local inexten- 
sibility of the chain (r'(s) 2 = 1) [55], the Hamiltonian 
reads 



H = 



1 



ds [ K r"( S ) 2 + /(s)r'( S ) 2 + 7 r±(s) ; 



(1) 



The polymer's overdamped motion thus follows from 
equating friction forces on the one hand and elastic and 
thermal forces on the other hand: 

CMM) = ^MH + {(«,t). (2) 

or(s, t) 

The slender shape of the polymer allows to decompose 
the friction matrix £ = d[(l — r' <E> r') + £ r' ® r'] into 
a parallel and transverse component [56], where the pa- 
rameter C 1/2 quantifies the friction anisotropy. 

In order to simplify notation, we shall set k = 1 
and £_l = 1 from now on, making forces a length -2 , 
times a length and the effective confinement strength 
7 a length -4 . Now, assuming the chain to be either 
strongly confined or sufficiently stiff, the transverse gra- 
dients = 0(e) <C 1 become very small (we will give a 
definition of e in terms of the system's parameters below) . 
Moreover, choosing the z-axis to coincide with the direc- 
tion of the polymer's extension allows us to effectively 
discriminate between transverse and longitudinal dis- 
placements and to write r(s) = (rj_(s),s — rii(s))*, where 
rn quantifies the amount of contour length stored in ther- 
mal undulations. For weakly bending polymers the in- 
extensibility constraint r' 2 = 1 implies rl — + o(e) 
and we may expand Eq. ^ in powers of e. Keeping 
terms up to order 0(e), we thus arrive at the following 
equations of motion, governing the chain's transverse and 
longitudinal dynamics: 



d t r ± = -i>»'- 7 rx + [fi' x ] , + tx 
(d t r l{ = (C - l)r'ld t rx - rff" - /' + 



(3a) 
£||(-3b) 
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Solving Eqs. ([3| by means of ordinary perturbation the- 
ory only works for large enough times, but runs into trou- 
ble in the limit of short times. We will thus present a 
refined perturbation theory in section [TV] 



III. SCALING PICTURE 



7 -V4 

£-1/2 

hit) 



Odijk's deflection length 
Euler buckling length 
(transverse) equilibration length 
boundary layer size 



TABLE I: Summary of characteristic length scales 



A. Length scales 

Much of the physics underlying the relaxation mecha- 
nisms of confined semiflexible chains can be understood 
by considering the interplay between the various char- 
acteristic length scales. Apart from the contour length 
L and the persistence length £ p , four additional length 
scales can be identified, which are crucial to the dynam- 
ics of such polymers (see table [I] for a summary). 

Two static length scales may readily be inferred 
from the externally applied force f and the confinement 
strength 7. A simple dimensional analysis reveals that, 
in our system of units, the length 7" 1 / 4 may be identi- 
fied with Odijk's deflection length [51] . whereas the scale 
f -1 / 2 , set by the externally applied force, is nothing else 
but the very length for which f represents the critical 
Euler load. We shall refer to this latter scale as "Euler 
buckling length" . 

Apart from these static scales, two characteristic dy- 
namical length scales have been identified to play a key 
role in the understanding of the polymer's relaxation dy- 
namics, the first of which can easily be deduced using the 
linearized version of the equations of motion ^ 



l/2-i 



Here we used 



/' = 



(4) 



(5) 



(to lowest order) from the expansion of the longitudinal 
part of ([3]). Denoting the characteristic length scale for 
transverse fluctuations by £ j_ , inspection of Eq. Q on a 
scaling level 



r 



(6) 



(where the externally applied force f sets the scale for the 
backbone tension /) immediately yields the scaling law 
for £ x 



£ x (t) 



i 1/4 , t <min{f" 2 ,7" 

vjt, r 2 «i«7 _1 - 



L } 



(7) 



As will be obvious from a more detailed treatment of the 
linearized equation of motion Q in section IV 



at time 



t the transverse fluctuations of the polymer are in equi- 
librium with their surroundings on length scales ~ £± (t) , 
which is why we refer to £± as (transverse) "equilibration 
length" . Note, in particular, that the second line of Eq. 



([7]) only applies in cases where the Euler buckling length 
is much smaller than Odijk's deflection length. 

This last statement can be re-expressed in more intu- 
itive terms. A simple dimensional analysis reveals that 
confinement induces a new force scale 7 1 / 2 which we will 
refer to as "effective confinement force" . The concept of 
identifying confinement as an effective force will be cen- 
tral in the further course of our discussion. In particu- 
lar, it allows to directly compare the strengths of applied 
forces and confinement. To this end, we introduce the 
"force scale separation parameter" 



7 



1/2 



(8) 



measuring the effective confinement force 7 1 / 2 in units of 
the externally applied force f. So for the second scaling 
regime in Eq. ([7]) to occur, the effective confinement force 
7 1 / 2 has to be small compared to the externally applied 
force f, viz. 

In the remainder of this contribution we shall mainly 
focus on the case x -C 1, where the chain's response to a 
large external perturbation is non-linear, discussing the 
case 1 > 1 of comparably small external perturbations 
along the way. The experimental accessibility of these 
limits will be investigated in section |VI[ 

According to the ordinary perturbation scheme em- 
ployed to arrive at the linearized equation of motion Q 
£j_ would be the only problem-specific dynamical length 
scale. As has been recognized before [3T)1 - |4"TI [57] , how- 
ever, ordinary perturbation theory (OPT) breaks down 
in the limit t — > 0, where longitudinal friction becomes 
important. As can be seen from Eq. ([5]), the effect of lon- 
gitudinal friction is of higher order in e and thus absent 
in leading order OPT. In order to account for the effect 
of longitudinal friction, at least on a heuristic level, the 
statement of a spatially constant tension, Eq. (fHl), has 
to be reconsidered. In the short time limit, longitudinal 
friction prevents the bulk of the polymer from respond- 
ing to sudden changes in ambient conditions and thus 
confines the longitudinal dynamics to boundary layers of 
dynamical size £\\(t) (cf. Fig. fll). Hence applicability of 
Eq. ([5| has to be restricted tolocal scales of size £±(t), 
with tension variations taking place at the much larger 
length scale £u(t). This separation of scales is the ba- 
sis for a more sophisticated multiple scale perturbation 
theory proposed in Refs. [43], [44] , which will be briefly 
reviewed in section Hvl 

Before we go on to discuss the scaling behavior of the 
chain's longitudinal response, we briefly introduce one 
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further central quantity in the discussion of the longitu- 
dinal dynamics of bio-polymers: the density of contour 
length stored in thermal undulations (p) = (r'n) s» |(r^}, 
where the last approximation holds for weakly bend- 
ing chains. Since we may apply equilibrium theory on 
scales ~ £x(t)t for the time being we are only inter- 
ested in the equilibrium value (p)o of this quantity, which 
we refer to as "stored length density" in what follows. 
To find an explicit expression for the spatial average of 
(p)o in the weakly bending limit, we expand the trans- 
verse displacements in simple sine-modes [55] rj_(s) = 
\Z^/L^2 n SL n sin(mrs/L). Noting that the stored length 
density reads (p} = (2L) _1 J2 n (nir / L) 2 (a.^) , where 
the fluctuation amplitude (a 2 ) can be determined by 
equipartition of energy, we find 



(p)o 



1 

H 



oo 



(T) 2 



v 



l(T) 4 + f(T) 2 +V 



(9) 



Depending on the magnitude of tension and confinement, 
this corresponds to the well known limits for the stored 
length density [5T1I55H5D]: 



(p)o 



6f,, 



. 2^271/2 



maxjf^ 1 / 2 } < L~ 2 
max{7 1 / 2 ,L" 2 } < f 
max{f, L~ 2 } < 7 1 / 2 . 



(10) 



Comparison of the stored length density for force- and 
confinement-dominated situations (i.e. the second and 
third line of Eq. ( |To| ) suggests to identify effective con- 
finement force and externally applied force (or more ex- 
actly 27 1 / 2 of). We will see in subsequent sections 
that this conclusion carries over to non-equilibrium situ- 
ations, where the dynamical scaling laws characterizing 
the relaxation processes for both force scenarios (pulling 
and release) can be identified with those of prestretched 
chains (cf. Rcf. [46 ). Most prominently we will re- 
encounter the effective confinement force in the scenario 
of free relaxation from confinement, where an initially 
tension free contour builds up a bulk-tension which is 
equal to 27 1 / 2 during relaxation in free space. More- 
over, Eq. ( 10 ) lends itself to define the small param- 



eter e quantifying the strength of transverse fluctua- 
tions, which was heuristically introduced in the last sec- 
tion. Recalling that (p) ~ (r^) ~ e, we find e £ 
{Lfp 1 , f _1 / 2 £p 1 , 7~ 1 / 4 £~ x }, depending on whether one is 
dealing with free, strongly prestretched or strongly con- 
fined chains, respectively. Note that the two latter defi- 
nitions of e allow for the use of the weakly bending rod 
approximation even for quite flexible polymers, such as 
dsDNA, provided the prestress or confinement is suffi- 
ciently strong. 



B. Longitudinal response 

Having identified all characteristic length scales, we go 
on to illustrate the physical mechanisms governing the re- 



laxation dynamics of harmonically confined semiflexible 
chains within the short time tension propagation regime 
t <C t| . Here the time scale t\ marks the time when the 
boundary layers eventually span over the whole contour, 
which is implicitly defined via ^(if,) = L. To this end 
we derive the scaling laws for the change in projected 
end-to-end distance 



(AiJ||)(*) = (JZ,|(*)-i2||(Q)>. 



(11) 



These results will be confirmed in the course of a more 
rigorous analytical treatment of the underlying equations 
of motion in section |IV| The arguments leading to the 
scaling picture are analogous in all three cases, whence 
we will confine ourselves to the discussion of the pulling 
scenario. 

In order to arrive at a scaling description for the chain's 
longitudinal response (Ai?||) we follow the reasoning in 
Ref. The boundary layers of size In, to which the 

longitudinal dynamics is confined, may notionally be sub- 
divided into much smaller segments of length £j_ <C £\\ , 
where the polymer's transverse fluctuations are in equi- 
librium with their respective surroundings. Each bound- 
ary layer is hence composed of £\\/£± segments and the 
scaling for (AiZn) may be inferred via the accumulated 
longitudinal response 



(Ai?||)(i) 



hit) 



(12) 



where the extension of each subsection 5 follows from 
equilibrium theory. The scaling for £n is determined by 
the observation that the force necessary to drag a poly- 
mer section of length £n through the solvent must not 
exceed the externally applied force f [31]. In mathemat- 
ical terms this condition translates to 



Chit) 



hit) s 
i±it) t 



f, 



(13) 



which may be inverted to yield the scaling of £\\ . 

Considering the pulling scenario, several asymptotic 
time regimes may be distinguished, depending on the 
magnitude of the equilibration length £± relative to the 
characteristic static length scales in the problem. For 
x -C 1 the Euler buckling length f -1 / 2 is much smaller 
than Odijk's deflection length 7 -1 / 4 , whence three dis- 
tinct asymptotic time regimes emerge within the tension 
propagation regime: 

1. For times t <C f~ 2 , the polymer is equilibrated 
on length scales much shorter than Euler buckling 
length (£j_ <C f~ 1//2 )- On local scales, the externally 
applied force f may thus be considered as small per- 
turbation and the chain relaxes freely (£± ~ i 1 / 4 ). 
The extension per segment S ~ (^l/^ P )f thus fol- 
lows from the linearized force-extension relation 
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2. For intermediate times f~ 2 « i « 7~ 1 / 2 f~\ the 
equilibration length clearly exceeds the Euler buck- 
ling length, but is still small compared to Odijk's 
deflection length (j^ 1 / 2 < t x < 7~ 1/4 )- Within 
this time regime the chain is yet not influenced by 
confinement on local scales. The externally applied 
force, however, may be regarded as large and the re- 
laxation behavior changes from free to forced relax- 
ation (£± ~ \fft). The small segments almost get 
completely stretched and we may estimate the ex- 
tension per segment 5 ~ £\j£ v by assuming all ini- 
tially stored length to be stretched out completely. 

3. Finally, once the equilibration length exceeds 
Odijk's deflection length at times 7~ 1 ^ 2 f~ 1 <C t <C 
t\, confinement suppresses long wavelength modes 
in the initial spectrum of the stored length density. 
Again, assuming this stored length to be entirely 
pulled out by virtue of the large force f we arrive 
at<5~( 7 - 1/4 /^)^- 

Having determined the local extension S within the 
various asymptotic time regimes, the scaling law for £\\ 
may readily be obtained by means of equation (13 1. 




' t l/8 
(fi)V4 

7 1 ' 



t « r 2 
r 2 < t < 
i 



(14) 



/ 7 f2 



< t < V 



Hence, using Eq 



find 




t « r 2 
r 2 < t < 



(15) 



/ 7 p 



It is intuitively clear that these scaling relations are 
identical to the pulling scenario for unconfined chains 
as long as £j_ <C 7 ~ 1//4 , i-e. as long as the chain has 
no perception of the confining channel walls on equili- 
brated length scales. In fact, the two time regimes for 
t <C 7~ 1 ^ 2 f~ 1 have already been discussed in the con- 
text of unconfined pulling in Refs. H3 HH US ES], 
where identical scaling relations for t\\ and (AiJii) have 
been discovered. In addition, analogous to the case of 
pulling chains subject to a prestretching force / pre [IB] , 
confinement introduces an extra asymptotic time regime 

(7 -1 / 2 f -1 < i < f|) whose respective scaling relations 
are identical to those of prestretched pulling upon identi- 
fying / pre o 7 1 / 2 . In this context confinement may thus 
be interpreted as effective prestretching force of magni- 
tude ~ 7 1 / 2 . 

Before we go on to a more rigorous treatment of the 
equations of motion (23), we briefly address the case 
x 3> 1, i.e. f -C 7 1 / 2 . Noting that, on a scaling level, the 
confined polymer may be regarded as being composed of 
L/7 -1 / 4 independent chain segments, the above state- 
ment x ^S> 1 means that these segments are very stiff 
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TABLE II: Summary of the scenario specific parameters oc- 
curring in Eqs. (25| and |28]). 



compared to the small externally applied force f, whence 
the whole chain responds linearly to the external per- 
turbation. Pulling and release may thus be regarded as 
mutually inverse scenarios. The Euler buckling length 
f -1 / 2 now exceeds Odijk's deflection length j~ 1 / 4 and 
only two time regimes may be distinguished during ten- 
sion propagation, depending on whether £j_ -C 7 -1 / 4 for 
times t < 7" 1 or £± > for times t > 7 _1 . The sit- 

uations for short times t <C 7" 1 as well as the respective 
scaling laws are identical to those for x <C 1. The scaling 
laws for t 3> 7" 1 may be obtained by noting that the local 
extension S ~ (£±/£ p ) 7 _3 ^ 4 f follows from the linearized 
force-extension relation in confinement [61] . Since the 
chain relaxes "tension- free" all the time (£±_ ~ i 1 ^ 4 ), we 
arrive at the following scaling laws for (A Jin) 




t < 7 ~ 
t > 7 ~ 



(16) 



Here the plus sign refers to pulling, the minus sign to 
release. In Eq. ( 16 1 we recover the well-known ~ t 7 / 8 
scaling law for short times and obtain a ~ i 1 / 2 scaling 
law for long times. Both scaling exponents have recently 
been deduced by Nam et al. |3S] for confined chains 
within linear response theory. 



IV. ASYMPTOTIC EQUATIONS OF MOTION 
A. Linearized dynamics 

Even in the limit of short times, where the action of 
longitudinal friction causes the breakdown of OPT, the 
linearized equations of motion Q remain applicable at 
local scales of size ~ £± (t) . A formal solution of Eq. Q is 
therefore the starting point for all subsequent discussions. 

Eq. @ may readily be solved by means of a Green's 
function method. Concentrating on bulk dynamics we set 
L — > 00 and determine the Green's function in Fourier 
space: 



where 



@(t - i) e -i 2 [i 2 (t-t)+F(t)-F(i)]- 1 (t-t) ^ q 7 -j 



F(t) 



dif(i) 



(18) 
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denotes the time integrated tension. Strictly speaking 
X± stated in Eq. (17) would have to be complemented 



where 



correspondingly, in order to account for problem specific 
boundary conditions. As discussed in Ref. [44] . however, 
the impact of boundary conditions is of minor importance 
for our present purposes, which is why we will dispense 
with any further discussion of this issue, referring the 
interested reader to Ref. [50] ■ Given the response func- 



tion x± m Eq. (17), mediating the influence of external 



and thermal forces, the mode spectrum of transverse dis- 
placements immediately follows by 



dtx±(q;t,t)$(q,i). 



(19) 



Later on we will need the coarse grained stored length 



density [57] (p) 
lows via Eq. (fl9 



(P>(*) = 



(r'D/2 + o(e), which immediately fol- 



dq 



o t4 



(p)o(q)x±(q;t,o) 



(20) 



+ 2q 2 J dt x ±(q;t,t)y 



where we used (£_i_(s, i)) = A/£ p S(s - s)5(t - t). 

Here the initial mode spectrum (p)o may be inferred from 
a continuum limit of Eq. (|9| . 



B. Multiple scale perturbation theory (MSPT) 



As briefly discussed in section |III[ the notion of a con- 
stant tension, as expressed by equation ([5]), has to be 
reconsidered on scales beyond the equilibration length 
£±(t). On the much larger scales ~ (\\(t) longitudinal 
friction causes the build-up of a non-trivial tension profile 
which in turn dominates the global relaxation dynamics 
of the chain. In order to understand these connections 
quantitatively, OPT from above has to be refined to take 
into account large scale tension variations. The key ob- 
servation, leading to an improved description of relax- 
ation dynamics suggested in Refs. [H], lies in the 
large separation between the typical scales characterizing 
transverse (£±) and longitudinal (£y) dynamics, holding 
true for weakly bending chains. This locally allows for 
the use of the linearized equations of motion. Having 
at hand the respective scaling laws governing £j_ and In 
(cf. section III), it is a straightforward matter to check 
that £x/(-\\ ;$ e *C 1 indeed holds in any of the cases 
considered in this work. As detailed in Refs. 03111], 
this observation may be exploited to conduct a multiple 
scale perturbation analysis, which eventually relates the 
curvature of the local tension profile to the change in the 
coarse grained stored length density: 



(Ap)(s,t) = (p)(s,t) - (p)(s,0). 



(22) 



In Eq. |2l| ) the stored length density (p) inherits its 
arc length dependence adiabatically from the (time inte- 
grated) backbone tension F. Using Eqs. ( |17[ ), (20) and 
(21 1 we thus arrive at the following closed relation all of 



our subsequent work is based on: 
d 2 s F(s,t) = 



dq 



1"*p 



(P)o(q) ( 



_ -2q 2 [q 2 t+F(s,t)]-2 7 t 



(23) 



2q 2 j rff e -2« 2 [« 2 (t-*)+-F(^t)--F(^t)]-27(t-t) 



In order to corroborate the scaling picture drawn in the 
previous section we now switch on to a rigorous analysis 
of Eq. (23). For the sake of convenience we introduce 



dimensionless variables which are defined by means of 
the following rescaling scheme: 



t 

q 

F(s,t) 




-1/4 



<r,T)-f 



(24a) 

(24b) 
(24c) 
(24d) 



Here f, denoting the externally applied force in pulling 
and release setups, needs to be replaced by 7 1 / 2 in the 
case of free relaxation from confinement. Applying this 
rescaling procedure renders Eq. ( 23 ) in the form 




_ -2q 2 [q 2 T+<f>(a,T)}-2bT 



q 4 + aq 2 + b 

rf - e -2 9 2 [ 9 2 (r-f)+0( C r,T)-0( C r,f)]-2b(r-f) 

(25)' 

where the scenario specific values of the parameters a and 
b are gathered in table |TT] All of our subsequent calcula- 
tions are based on (intermediate) asymptotic expansions 
of Eq. ( 25 ) , assuming the time integrated backbone ten- 
sion 4> to attain the following one parameter scaling form 



(26) 



with scaling exponents a and 77 depending on the par- 
ticular asymptotic time regime under consideration. The 
validity of this scaling assumption will be confirmed in 
the following subsections, where the values of a and -q 
are explicitly calculated. We also note in advance that 
the applicability of our subsequently developed approxi- 
mation schemes relies on the ad hoc condition 



d 2 s F(s,t) = -C((A P )(M), 



(21) 



a>-l/2, 



(27) 
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which will indeed be justified with the benefit of hind- 
sight. We complete our mathematical description of the 
system by stating the following scenario specific bound- 
ary conditions for semi-infinite filaments, supplementing 
Eq. 



0(0, r) 



co 



lim D 4>{a, t) = c c 



(28a) 
(28b) 



Here the scenario specific constants cq and Coo and the 
operator D are given in table [XT] The assumption of deal- 
ing with semi-infinite contours amounts to concentrating 
on the dynamics at one of the polymer's ends only. This 
is possible since, during the tension propagation regime, 
the dynamics entirely occurs within boundary layers of 
size £|| (t) <§C L which are separated by a large bulk sec- 
tion and thus are independent from each other. For later 
times, of course (i.e. when £\\(t) is comparable to the 
contour length L), this assumption has to be reconsid- 
ered. 



since the term multiplying this exponential vanishes for 
wave numbers q > r -1 / 4 3> 1, where this constitutes a 
bad approximation. We are thus left with 



io H 1 



(30) 



v 



dr' cj)(a,T r )e- 2q (r ~ 7 



Finally, mapping Eq. ( 30 ) to Laplace space and using 



x <C 1, we eventually arrive at the following linearized 
asymptotic differential equation 



(31) 



which is to be solved subject to the boundary conditions 
(cf. Eqs. pi)) 



0(0, z) = z~ 
lim 4>(<j 1 z) = 0. 



(32a) 
(32b) 



C. Pulling 



Adjusting Eq. (25) to the pulling scenario (i.e. setting 
a = 0, b = x 2 <C 1) we seek possible asymptotic expan- 
sions for the various time regimes identified in section |Hl| 
Here the boundary conditions, stated in Eq. (28), reflect 



the action of a constant pulling force f at the ends, and 
a vanishing backbone tension in the bulk of the chain. 



1. Short times 

We start our discussions with asymptotically short 
times r <C 1. Invoking our general scaling assumption 
Eq. (26) and assuming a > —1/2, we conclude that 
the bending modes, which are of order 0(q 4 r) domi- 
nate over tension modes 0(q 2 T a+1 ) in the exponentials. 
This allows us expand the exponentials up to first or- 
der in the tension modes. Since, moreover, we cur- 
rently concentrate on cases with x <C 1, we may set 



-2x z t 



1 



a x (t— f) 



arriving at 




The first factor in line one of this equation is of order 
0(q 4 r) -C 1 for q < t -1 / 4 and bounded for q — » oo, 
whence the entire first line may be neglected. Moreover, 
the exponential occurring in line two may be set to unity, 



In Eq. (31) <fi(a,z) = J^dre zt 4>(<j,t) denotes the 
Laplace transform of 0(er, r). One easily verifies, that 



/ \ -2 -2- 3/s z 1/8 a 



(33) 



is the solution to the boundary value problem stated in 
Eqs. (31) and (32). The real space solution, following 



from Eq. ( 33 ) by means of a Laplace back-transform 



£+MO dz^ 

27T* ./e— ioo z2 



-2- 3 / s zV 8 <r+zr 



e'+ioo 
e' —zoo 

a 



dz _ 2 - 3 /8 z l/8_iZ_. 

— e t 1 ' 8 



(34) 



confirms our scaling assumption Eq. ( 26 ) with a = and 



7] = 1/8. Here we substituted zt — > z in the second line. 
Identifying 



T'l 



we immediately read off 




V jl/8 



(35) 



(36) 



corroborating our previous scaling analysis. As already 
discussed on a scaling level in section |III[ the short time 
asymptote of Eq. (25) (with a = and b = x 2 ) is iden- 
tical to the scenario of unconfincd pulling. We thus skip 
the explicit calculation of (f>(a, r) in Eq. (34), referring 
the interested reader to Rcf . [45 , where a more detailed 
discussion of this function, including a summation for- 
mula for <f>, is given. 



2. Long times 



below), we thus arrive at 



In order to discuss the situation for long times f 2 i 
t 1, we split the right hand side of Eq. (25) into 



a deterministic (D) and a stochastic (S) contribution. 
Adjusting the parameters a and b to the case of pulling, 
we write 



2 d _ -2q 2 [q 2 T+^(a,r)}-2x 2 



D 



S = - 



dq 9 2 (l - e 



q 4 + x 2 



" ^ 2q 2 I dfx 



(37a) 
(37b) 



-2q 2 \q 2 (T-f)+c/>(o-,T)-ct>(<T,f)]-2x 2 (T-f) 



Starting with the stochastic contribution S, we first ap- 
ply a quasi-static approximation, assuming the chain to 
be equilibrated under the current local backbone tension 
[42| . This technically amounts to an expansion of the in- 
tegrated tension <^>(<7, f) in the exponential about f = r, 
i.e. approximating 4>{t) — 4>{f) ~ <9 r 0(r)( T — t). Invok- 
ing our general scaling assumptions (26) and p7|) this 



strategy is rigorously justifiable [45 . As a result, we are 
left with 



dq 1 — e 



-2 \q 2 d T 4>(a,T)+x 2 T 



<7 2 + d T <j){a,T) 



(38) 



Here we omitted the term x 2 <C 1 in the denominator 
and invoked dominance of tension modes to neglect the 
bending modes in the exponential. Eq. ( 38 ) may even 



be simplified further. To this end we note that the func- 
tion (q 2 + drcj))^ 1 decays over a distance Aq = 0(r Q / 2 ), 
whereas the function (1 — e~ 2 ['"J r ) approaches unity 
within a much smaller interval Sq — 0(r~( a+1 '/ 2 ) <C Aq. 
We may therefore set the latter function to unity and find 



5' 



dq 



1 



-1 



7T q 2 + d T (j)(a, t) 2y/d T (f){a,T)' 



(39) 



In contrast to the stochastic contribution S, whose 
asymptotic form (39) follows for all times r 3> 1 from 



a quasi-static approximation, the deterministic contribu- 
tion D possesses two distinct asymptotes in the long-time 
regime. To see this, we make use of dominance of tension 
modes for r ^ 1 to write 



D 



dq1_ 

„ 7T 



1 _ e -2q 2 cj>(<j,T)-2x 2 T 



q 4 + x 2 



(40) 



Now the factor (1 — e~ 2<? ^~ 2x T ) may be set to unity for 
all wave numbers q > q* = t~(" +1 )/ 2 , whereas the ratio- 
nal function it multiplies reaches its maximum at q rn — 
x 1 ' 2 . We thus have to keep the factor (1 - er 2 i 2 <t>-2x 2 r- ) 

if q m < q,.Le. if t < 7 -l/(2o+2)f-(2a+D/(a+l) ] but may 

omit it once the opposite is true. Anticipating a = (see 




r 2 < t < 



\Af 2 (41) 



q 2 -\-x z 



hf 2 



< t < t 



L • 



where we have to keep the small term x 2 in the second 
line in order to avoid artificial divergencies for q — > 0. 
Note, moreover, that x 2 t — jt -C 1 for t <C (7f 2 ) -1 ^ 2 
and 7 f 2 , enabling us to approximate e~ 2x T ss 1 in 



the first line. Performing the integrals in Eq. (41) we 
obtain 




r 2 < t < 



V7f 2 L 



(42) 



As is obvious from physical considerations, the actual 
tension along the contour is comparable to the externally 
applied force, whence d T cj) = 0(1). This implies S — 
0(1) and the deterministic term (D = 0( T ( Q+1 )/ 2 ) > 1 
and D = 0(x~ 1 ^ 2 ) ^ 1 respectively) clearly dominates 
the dynamics in both intermediate asymptotic regimes. 
Thus, neglecting stochastic contributions and setting 



d 2 M°,T)~D, 



(43) 



we eventually arrive at the sought-after asymptotic equa- 
tions of motion. 

Here the intermediate asymptote D oc c/) 1 ^ 2 reproduces 
the "taut string approximation" , first used in Ref. |39j on 
a heuristic level in order to discuss tension propagation 
for large longitudinal pulling forces in free space. In the 
second long time asymptote, the deterministic contribu- 
tion D cx x~ x l 2 is nothing but the initial stored length 
density (in rescaled units) for polymers equilibrated un- 
der zero tension in cylindrical confinements of constant 
strength 7. Hence, the corresponding equation of motion 
might have been readily inferred from our scaling picture 
of section |TTT| by assuming the initially stored length to be 
pulled out completely by virtue of a strong pulling force 
f- 

As anticipated in section |III| the equation of motion 
governing the tension dynamics in the regime f~ 2 <C t <C 
(7f 2 ) -1 / 2 , where the chain yet has locally no perception 
of the confining channel walls, is just the same as the 
one found in Ref. |45j for unconfined pulling and asymp- 
totically long times. There the corresponding differential 
equation is discussed in detail. We simply state the so- 
lution: 



for) 



1 

72^ 



1/4 



£-1 



(44) 



again corroborating our scaling assumption with a 
and 77 = 1/4 and hence 




(45) 
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The second asymptote in Eqs. (42) and (43) may be 



written down in a parameter free form by means of a 
redefinition of the dimensionless arc length variable a. 
Setting 



(46) 



and invoking our scaling assumption ( 26 ) with 77 = 1/2 
and a — we find 



dU(0 



1 



2^2' 



(47) 



implying a parabolic ^-dependence for </>. Noting that 
the actual backbone tension is given by / = f ip, where 
tp = (f> — £d£<fr/2, we see that Eq. (47) predicts a lin- 
early decreasing tension profile. Obviously, this result 
is meaningful only within the boundary sections of the 
chain, since a linearly decreasing tension profile cannot 
smoothly match the boundary condition in the bulk Eq. 
( 28b ). Of course, our above reasoning leading to Eq. 
(43) was based on the assumption (p = d T (f> = 0(1), 



which certainly breaks down far away from the bound- 
aries, where the local backbone tension <gC 1 is close to 
zero and stochastic effects dominate the physics. A more 
complete description thus has to include the influence of 
thermal motion in order to account for the smooth tran- 
sition between a linearly decreasing tension within the 
boundaries and a constant (zero) tension in the bulk |46j . 
Since the details of this transition, however, are dispens- 
able in the context of longitudinal dynamics, we will not 
enter this discussion and confine ourselves to the simple 
picture suggested by Eq. ( 47 1 , predicting a piecewise lin- 



ear tension profile. A comparison of this piecewise linear 
estimate to a full numerical solution of Eq. ( 25 ) , reveals 
that this is indeed a feasible approximation. Using the 
scaling exponent 77 = 1/2 and recalling the redefinition 
of a in Eq. ( 46 ) we infer 



/u(t)*y|7 1/8 (fe) 1/a , (48) 

as expected from our scaling picture. 



OPT regime 



So far we concentrated on short times t <C t L , where 
the bulk of the polymer is not influenced by the externally 
applied tension f, and where we were able to restrict our 
discussions to formally semi-infinite chains. In contrast, 
for times t ^> v\ tension propagated through the whole 
contour and we have to account for the finite length of 
the chain in general. We estimate the crossover time t\ 
by setting ^(tfj = L. Using Eq. @ we find 



T 
v 



h C l v 7 1/4 f 



(49) 



We shall now show, that this time scale is indeed identical 
to the time i*, marking the crossover to applicability of 
OPT. To estimate note that OPT predicts a flat ten- 
sion profile /opt(s) = f, which in our rescaled version of 
the time integrated tension amounts to 4>opt{o~, t) = r. 
In reality the tension spatially varies along the con- 
tour and we may write 4 > { <J ^ t) — t + S(f>(a, t) for times 
t = 0(t*) @B]. Redefining a = s/L and using Eq. ((47}, 



which is still valid for times t > t L , we find the following 
boundary value problem governing 5 (jr. 



I? f 



l P 7 1/4 



5<j>(a,r) 



<T=0,1 



(50a) 
(50b) 



The solution to this problem may readily be inferred 
by straightforward integration. We immediately read off 

5(j) = O y^r). In order for OPT to be applicable 

we have to require 5<j) <C </>opt = t. From this we deduce 
the time scale signaling the crossover to OPT 



L 2 



l v 7 V*f 



(51) 



which turns out to be identical to t L stated in Eq. ( 49 1 



Hence we are in a position to invoke OPT in order to 
determine the longitudinal dynamics in the time regime 
succeeding tension propagation. 



D. Release 

It turns out that confinement effects remain invisible 
to the relaxing contour during the entire tension prop- 
agation regime, due to a large bulk tension ~ f, which 
suppresses thermal undulations to amplitudes far below 
the channel diameter. Within this tension propagation 
regime, we will therefore confine ourselves to a brief dis- 
cussion of the approximations necessary to map the equa- 



tion of motion (25) (with a = 1 and b = x ) to the re- 



spective asymptotic differential equations, which indeed 
turn out to be identical to ones discussed in Ref. [45] in 
the context of unconfined release. The actual presence 
of confining channel walls becomes visible not before the 
tension has essentially relaxed, resulting in a correspond- 
ing shift in the time scale i*, which marks the crossover 
to OPT (cf. section [TVTT2|. 



1. Tension propagation 

The universal short time regimes t«1 may be treated 
in an analogous fashion as above (cf . section IV C 1 1 . In- 



voking dominance of bending modes and linearizing Eq. 
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(52) 



Noting that only wave numbers q ~ r^ 1 / 4 3> 1 contribute 
in the first line, we may set q 2 (q 2 + 1) ~ q 4 . Moreover, 
the exponential occurring in line two may be set to unity 
again. Mapping the resulting expression to Laplace space 
and performing the remaining g-integrals we are left with 



< 2 J(a, z) ~ (2- 3 /V/4 - l) 0>, Z ) - 2- 3 ^z-y\ (53) 



subject to the boundary conditions (cf Eq. (28)) 

0(0, z) = (54a) 
lim 4>{a,z) = z~ 2 . (54b) 

a— ► oo 

For an explicit treatment of this particular boundary 
value problem in the short time limit z » 1 we refer 
the interested reader to Ref. [15]. There, again, our pre- 
viously proposed scaling form Eq. ( 26 1 is confirmed with 



a = and rj = 1/8 giving rise to the universal short time 
scaling law 




(55) 



in accordance with our scaling picture. 

The asymptote for long times r> 1 may be inferred 
along exactly the same line of reasoning as outlined for 
the pulling scenario, giving (x <C 1) 




1 



-2q 2 d T 4>-T 



d T <t> 



(56) 



Moreover, using d T <f> = 0(1), the rational functions (q 2 + 
l)" 1 and (q 2 + 9 T 0) _1 decay over distances Aq — 0(1). 
The exponentials, in contrast, decay over much smaller 
distances Sq = 0(T~ (a+1 *>/ 2 ) -C 1. It is thus legitimate 
to neglect the exponentials altogether, to arrive at the 
following asymptotic differential equation 



(57) 



governing the system's tension dynamics in the long time 
regime. Once again, choosing a — and r\ = 1/2, this 
equations may be solved by means of a one parameter 
scaling function of type Eq. ( 26 I , giving in particular 




(58) 



2. Homogeneous tension relaxation and crossover to OPT 

Unlike in the case of pulling, the time scales t| and t* 
are separated by a regime of homogeneous tension relax- 
ation 45 in the release scenario. In particular, setting 

£u(t[) = L we find 



f \\ _{ L ^c-3/2 



(59) 



which is identical to the respective crossover scale ob- 
served in unconfined release 45 . In contrast to its un- 
confined counterpart, however, the crossover to applica- 
bility of OPT, where the backbone tension has essentially 
relaxed to its new equilibrium value f(s,t) = 0, occurs 
significantly earlier since confinement supersedes relax- 
ation of modes q < 7 1 / 4 . To estimate the time scale t* 
marking the transition to OPT we calculate the change 
in stored length predicted by OPT and compare the cor- 
responding friction force (L 2 (Ap)oPT/t to the charac- 
teristic force of the currently relaxing mode. Noting 
that the right hand side of Eq. ( [25| i s proportional to the 
change in stored length (cf. Eq. (|2ip), we obtain (setting 
= 0) 



(Ap> 



-1 



OPT 



o 



dq 



1 



q A + x 2 



-1/4 



(60) 



which could have been anticipated by means of equilib- 
rium theory. In deriving this result we used 1 < 1 and 

omitted the factors ^1 — e~ 2 ( < ? 4 + :r2 ) r ^ ; which approach 

unity within a fairly small interval 8q = 0(r -1 / 4 ) <C 1 
and thus do not affect the integrals. Comparison of fric- 
tion and characteristic force (J 2 ~ t -1 / 2 (chain relaxes 
freely, i. e. under vanishing tension) thus provides us with 
the crossover time scale 



r 4 



(61) 



which — as anticipated — is much smaller than the respec- 
tive time scale for unconfined chains, provided the con- 
finement is sufficiently strong (cf. section [V]) . Given that 
x <C 1 and f > f c the time £* is indeed large compared 
to the crossover scale i'l 



1/2 



X' 1 > 1. 



Here 



i 2 



C 2 £ 4 



(62) 



(63) 



denotes a critical scale for the externally applied force, 
below which the chain reacts linearly to the external 
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perturbation, never entering any non-linear long-time 
regimes. Focusing on the more interesting cases of non- 
linear response we tacitly assume the inequality f > f c to 
be fulfilled. Under these conditions a regime of homoge- 
neous tension relaxation emerges for times t \ <C t <C , 
whose dynamics, governed by the same equation of mo- 
tion (57 1, is identical to the one within the respective 



unconfined scenario, where the tension is shown to decay 

~ r 2 / 3 051. 



E. Free relaxation from confinement (FRC) 



transform the resulting expression. Performing the re- 
maining g-integrals and focusing on short times, i.e. 
z» 1, we arrive at 



8 2 J(a, z) ~ 2- 3 /V/ 4 0(a, z) + 2- 1 l 2 z~\ 



(65) 



which is subject to the boundary conditions (cf. Eq. 



0(0,*) = 
lim d fT cf)(a, z) = 0. 



(66a) 
(66b) 



As a third scenario we shall investigate the intermedi- 
ate asymptotics in the case of an initially confined con- 
tour relaxing in free space. While strongly increasing 
the confinement strength eventually leads to creation of 
hairpin structures, analogous to the case of compressing 
a polymer .62] , and thus forbids the use of the presently 
applied theory for late times, a sudden drop in confine- 
ment strength is devoid of such peculiarities and may be 
described successfully by means of Eq. ( 25 1 . Here we 



restrict ourselves to the discussion of a sudden change in 
confinement strength from some finite value 7 to zero. 
This scenario may be regarded as another typical exam- 
ple of initially straight contours relaxing in free space 
[50| . We will show that this problem is indeed very sim- 
ilar to the unconfined release scenario [151 [50], safe for 
a transient short time regime, where tension builds up 
along the contour and the longitudinal dynamics (AiZii ) 
will be shown to follow a new superlinear scaling law. In 
addition, the prefactor for (AiZii) in the long time ten- 
sion propagation regime turns out to be different, due to 
a somewhat different condition in the bulk. 

As a key result of the present subsection we will deter- 
mine the scaling exponent quantifying the growth of the 
backbone tension at short times. A quantitative discus- 
sion of the longitudinal response (Ai?||) will be given in 
section |Vj 



This linear second order boundary value problem is 
solved straightforwardly by standard methods yielding 



( ( T,z) = 2 1 /4 z -9/4^ 1 _ 



e r"'* 1 "*) . (67) 



From this, the sought-after solution <p(cr, r) is found by 
means of an inverse Laplace transform. Using the linear- 
ity of the inverse Laplace transform, as well as 



c- 1 



y -9/4+n/8 



= T 5/4 



-n/8 



r(S-f) 



(68) 



where n ^ 18 + 8A; (k = 1, 2, . . . ) and denotes the 

inverse Laplace transform, we may write down <j)(cr, r) in 
the form of a summation formula: 



2 l/4 £ -l L-9/4 A 

XJ ( 1 \n+lo-3n/8_n 

2^4 tU f £-1 

n=l 

9 i/4^5/4 ^— J t n 



-9/4+n/8 



1 -T(f-t) 



(69) 



1. Short times 



For short times (r <C 1) we again linearize Eq. (25) 
(with a — and 6=1): 



dq J e- 2 ^ T - 1 
V] <Z 2 (<? 4 + 1) 

+ 20(a,r) \l + e^ T (^rT~l)] «>4i 

4q 4 [ dr' 4>{(T,T l )e 
Jo 



Applying the same reasoning as before, we set the expo- 
nential in line two equal to one and subsequently Laplace 




Here we expanded the exponential in the first line and 
exclude the values n= 18 + 8fc (k = 1,2,...) in the sum- 
mation range. This result corroborates our general scal- 
ing assumption Eq. (26) with a = 1/4 and r\ = 1/8. The 
latter exponent again yields the known universal short 
time scaling law £11 ~ i 1 / 8 . The asymptotic form for the 
rescaled bulk tension ipoo = foo/f may readily be inferred 
from the bulk limit of Eq. ( |67| ). Noting that Cp = z<f> we 
conclude (poo = linv^oo zip = 2 1 / 4 z~ 5 / 4 . Thence, using 
Eq. (f68b we find 



2 1 ' 4 t 



^coM^^yT^ 4 . (70) 

The backbone tension thus grows ~ t 1 / 4 for short times. 
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2. Long times, homogeneous tension relaxation and OPT 



F. Pulling and release for small forces 



Following our approximation schemes explained for the 
two force cases Eq. ( 25 ) reduces to 



d 2 Ma lT ) 



dq 



1 



o * U 4 + l 



- d T 4>{(7, T) 



(71) 



for long times x> 1. Taking the bulk limit of Eq. (71) 



and using Eq. (28b) supplies us with the rescaled bulk 
tension. We fine 

=i~vibr^= 2 ' (72) 

where we used lim cr _ > . 00 9 T 0(cr, r) = ^oo(r). The bulk 
tension thus saturates at a constant value 



- 9^/2 



fee =2 7 



(73) 



for asymptotically long times during tension propaga- 
tion. Note that the asymptotic equation of motion for the 
rescaled tension ip — d T (j), entailed by Eq. (71 ), is identi- 



cal to the respective equation for unconfined release, safe 
for a shift in boundary conditions: in the present case 
the rescaled tension tp approaches the bulk value ifoo — 2 
rather than ip^ = 1. This in turn gives rise to a new pref- 
actor in the scaling law for (AiZn), which we determined 
by means of a shooting method (cf. table |III| and section 
V|). Again the scaling for £\\(t) is obtained by mapping 
q. ( 71 ) to an ordinary differential equation by means 



of our general scaling assumption ( 26 1 and identifying 
£ = s/iii (i). Choosing a — and 77 = 1/2 we find 




P 7 3/8 t l/2 j 



(74) 



which is the well known result from unconfined release 
[4"5] (identifying 7 1 / 2 -n- f). 

As soon as the boundary layers span over the whole 
contour and we have to take into account the finiteness 
of the chain. The mathematics describing the polymer 
dynamics now becomes strictly identical to the uncon- 
fined release scenario, enabling us to carry over the cor- 
responding results within the regimes of homogeneous 
tension relaxation and OPT. This in particular gives the 
crossover time scales: 



4 = cfV 3/4 

** — p4 ' 



(75a) 
(75b) 



which are again separated provided the effective confine- 
ment strength is sufficiently strong. 



For small forces f <C 7 1 / 2 , i-e. x 3> 1, pulling 
and release become mutually inverse scenarios. Since 
the respective short time asymptotes remain unchanged 
we restrict our discussions to asymptotically long times 
t > 7 _1 during tension propagation. The equation of 
motion |23l in the case of pulling reads 



d 2 s F(s,t) 



C r°°dq f q 



tpJo t U 4 + 7 q 4 + q 2 d t F(s,t)+~f 
-1/4 C f°°dqi q* ql 



(76) 

where we invoked a quasi-static approximation (which is 
justified at any rate, since jt 3> 1) and where we rescaled 
wave numbers q — > q^f 1 ^ 4 - Note that, in order to describe 
the release case, we would simply have to add the term 
q 2 x _1 in the denominator of the first term. However, 
since this term contributes only for wave numbers q > 
x 1 ! 2 where the complete first term is close to zero, it is 
justified to neglect this contribution altogether, whence 
pulling and release are both described by means of the 



same Eq. @. Noting that 7" 1 / 2 9 t F(s, t) = Oix- 1 ) < 



1, we may expand the integrand in Eq. (76) and find 



the following diffusion equation for the time integrated 
tension 



d t F(s,t) ~8V2-^j 



^ s 2 F( S ,t) 



with diffusion constant 



A oc 



7 3/ % 
C 



D~ 



(77) 



(78) 



in our units proportional to the inverse cross-sectional 
area of the channel (D denoting the channel diameter). 
Here we used Odijk's scaling law j~ 3 / 4 ~ L\ ~ D 2 l p 
|51j . From this we immediately deduce the scaling for 
the boundary layer size L\ [55] : 



^ J|7 3/ V /2 ocyAo7. 



(79) 



Once again, Eq. (77) may be solved by means of a one 



parameter scaling function. Defining 

_ f 

a = s 

V8V2Ao 



r = ft 



(80a) 
(80b) 



and setting F(s,t) = f 1 (f>(cr, r) we obtain 

d T ^(a,r)=d 2 a( f > (a,r). (81) 
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Or, using our general scaling assumption Eq. ( 26 ) with 
a = and r) = 1/2 we get 



(82) 



Hence, invoking the respective boundary conditions Eqs. 
(28), we find the following one parameter solution for 



pulling: 



0(C) - 1 - erf 



(83) 



scenario. Thus knowledge of the previously determined 
dynamical exponents a and r/ immediately supplies us 
with the scaling laws governing the longitudinal response 
(Ai?,||) within the respective asymptotic time regimes. 

The only unknown quantity c%<^(£ = 0) can be deter- 
mined from the corresponding equations of motion. We 
determined this numerical constant explicitly for the var- 
ious scenarios and time regimes discussed in this work, 
either analytically on the basis of the respective asymp- 
totic equation of motion or numerically by means of a 
shooting method. 



and for release: 



HO = erf 



(84) 



Here erf(£) denotes the error function. 

For sufficiently weak forces, the relaxation dynamics 
directly crosses over to the OPT regime in both cases. 
To estimate the corresponding crossover time it is thus 
sufficient to set £\\(t*) = L giving 



(85) 



T 2 



V. LONGITUDINAL RESPONSE 



Tension propagation and homogeneous tension 
relaxation 



Up to now our quantitative discussions have focused on 
the analysis of Eq. (25) and the derivation of respective 



scaling solutions for the tension profiles. These results set 
the stage for an analytical discussion of the longitudinal 
response of the chain, represented by the change in the 
projected end-to-end distance (AR\\)(t). This quantity is 
formally expressed in terms of the stored length density: 

(AR l{ )(t)= [ ds(rUs,0)-rUs,t)) 

rL (86) 
ds{Ap)(s,t) + o(e). 



Omitting terms o(e) and using Eq. (21) we may thus 
relate longitudinal response and (time integrated) tension 
profile 



<AJ2[|>(t) = i f dsd 2 F{s,t) = -^d s F(0,t) 
(Jo C 



-2dtf{Z = 0) 



c5/4+2(a 



-v) 



.a+1— 7) 



(87) 



where we used the symmetry of the boundary conditions 
and explicitly inserted our scaling assumption Eq. ( 26 1 . 

7 1 / 2 in the FRC 



Recall that we have to replace f 



B. Ordinary perturbation theory (OPT) 

Within the OPT regime, where the tension profiles are 
flat, the longitudinal response is calculated directly from 
change in stored length density. In the following we give a 
brief derivation of the respective scaling laws. To this end 
we calculate the change in stored length density (Ap)(t), 
which is proportional the right hand side of Eq. ( |23[ ). 
Within OPT this expression considerably simplifies due 
to the spatially constant tension profiles. Restoring orig- 
inal units and starting with the pulling scenario we find 



(Ap)op T (t) - 
dq 

dq 



1-e 



-2q 2 ft 



q 4 + 7 <7 4 + <7 2 f + 7 

9 

T 



l 7- 1 / 4 



Here we invoked dominance of tension modes and even- 
tually neglected the exponential in the second line, which 
decays on a scale of the order 0((f t) -1 / 2 ) which is much 
smaller than the scale characteristic of the rational func- 
tions (C(7 1/4 )) for times t > 7 -1 / 2 f _1 . Thence, noting 
that x <C 1, we find the change in stored length density 
to be identical to the (negative) stored length density of 
a polymer equilibrated in a constant confinement of effec- 
tive strength 7, corroborating our scaling picture drawn 
Finally, using (ARu) = 



in section 
at the sea 



III 



(Ap) L we arrive 

mg law given in table |III| 
For the FRC scenario we have to keep the exponential 
in the stochastic term in order to avoid artificial diver- 
gencies for q — > 0. We find 



(Ap)op T (i) 



93/4 

where we used t 1 / 4 7 -1 / 4 . 



dq J q 2 1 - er 2 i H 




(89) 
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Time Regime 


Pulling 


Release 


Free Relaxation From Confinement 


short times 


r 2 5/8 /pc 15 /g\l £ 7 / 8 


- f2 5/8 /re 15 /8ii (f/\flin) t 7/a 


-[2 7/8 /rfl7/8)l (V /4 / \ Kin) t 9/8 


intermediate times 


To //^n \l/4l/r3/4/ //^ \ ,3/4 

[8/(72tt) 7 ](f ' /\jC,t p )t ' 






long times 


2 3/4( f l/2 7 -l/8/y^) 




-2.946( 7 1/8 /y / ^) t 1/2 


t{ < t < u 










2-3/2^-1/4^) 


_ 2 -3/ 2(7 -i/4 L/ ^ ) 


-^/rii/mi/e^t 1 / 4 



TABLE III: Response of the projected end-to-end distance (AR^)(t) for large forces (x <C 1). The numerical prefactor for the 
scenario of free relaxation from confinement (long times) was newly determined by means of a shooting method. The prefactor 
for pulling (long times) was determined analytically on the basis of Eq. (471 by requiring the time integrated tension to be 
continuously differentiable at the transition to the bulk. 



Finally, the change in stored length density for 
release scenario was already calculated in Eq. ( 60 ) . 



the 



We conclude this section by briefly discussing the case 
of large force scale separation parameter x 3> 1 in the 
context of pulling and release. Note that within OPT, 
pulling and release are mutually "reverse" scenarios (cf. 
table III ) . Using equation ( 88 1 and rescaling wave num- 
bers (q — > we find 



(Ap)(t) ~ ± 7 - 1/4 



dq 



ir£ p [ql + x-hf + l q 4 + l 
dq q 4 



o nip ( q 4 + !)- 



+ 0{x- 2 ) 



2 -7/2 f 



tp 7 3/4 ' 



(90) 



where we expanded the first line in x^ 1 <C 1. In Eq. (90) 
the upper sign refers to pulling, the lower to release. 

Our results are summarized in tables \UJ\ and [PV] Here, 
as already discussed in the previous section, the force sce- 
narios are identical to their unconfined counterparts for 
sufficiently short times [45], whereas the dynamical scal- 
ing relations are analogous to the cases of prestretched 
chains for late times |46j , consistent with our earlier con- 
clusion to interpret confinement as effective prestretching 
force. Unlike the respective free space scenarios, how- 
ever, the force scenarios in confinement saturate at a 
constant value for (Ailii) within the OPT regime t > J,, 



Time Regime 


Pulling / Release (x 3> 1) 


t < 7- 1 


±[2 5 / 8 /r(i5/8)](f/y / c^)t 7 / 8 


7" 1 « t < 4 




t\=U<SZ.t 


±(2- 7 / 2 L/£ p )f/7 3/4 



TABLE IV: Response of the projected end-to-end distance 
(A_R||)(t) for the pulling and release scenario in strong con- 
finement for weak forces (i> 1). The corresponding scaling 
laws only differ by a sign and are linear in the externally ap- 
plied force f. 



which is due to the suppression of long wavelength modes 
by confinement. The FRC scenario bears close resem- 
blance to the scenario of unconfined release for late times 
t ^> 7 _1 , safe for a shift in the prefactor of (Ai?y) which 
is due to a somewhat different boundary condition in the 
bulk (cf. previous section). In particular the univer- 
sal (Ai?y) ~ i 1 / 3 scaling, characteristic of freely relaxing 
"initially straight" contours (cf. Ref. [50]), is recovered. 
Most remarkably, the FRC scenario exhibits a superlin- 
ear (A-Rn) ~ t 9 ^ 8 scaling for short times, which may be 
traced back to a steadily growing tension in the bulk. 



C. Weak confinement 

Our hitherto discussions for large external forces (x <C 
1) revealed that the two force scenarios are governed by 
the dynamical laws of the respective unconfined problems 
for sufficiently short times. In the limit 7 — >■ 0, of course, 
all the confinement induces differences vanish and we be 
able recover the results of Ref. [45] . We will now show 
that there exists a finite critical value of the confinement 
strength 7, below which confinement effects are irrelevant 
for the longitudinal dynamics. 



1. Pulling 

The pulling scenario starts to deviate from its un- 
confined counterpart as soon as the equilibration length 
£j_(t) grows beyond Odijk's deflection length 7 1 / 4 at 
times t > 7 ~ 1 / 2 f. Subsequently the chain's dynamics 
is governed by scaling laws which closely resemble those 
of prestretched polymers |46j and have no analog for un- 
confined molecules without prestress. In contrast to un- 
confined chains the longitudinal dynamics saturates at a 
constant value of the projected end-to-end distance for 
times beyond i*, which is thus identical to the longest 
relaxation time. 

In order for these observations to be true, however, we 
tacitly assumed the following ordering in time to hold: 
-1/2^-1 ^ ^ where is given in Eq. (51 1. 



r\<r 

While the first inequality is true by virtue of a large force 
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separation i« 1, the second inequality imposes a con- 
straint on the confinement strength 7. More precisely we 
have to stipulate 



7 > 7c = f c 



(91) 



The critical force f c (Eq. 



( 63 ) ) thus translates into a 
below which the third 



critical confinement strength 7 C 
asymptotic time regime within tension propagation van- 
ishes. Violation of Eq. (91) has even further conse- 
quences. First of all, since the scaling law lu ~ (ft) 1 / 4 
(rather than In ~ (ft) 1 / 2 ) governs the growth of l\\ at the 
crossover between tension propagation and OPT regimes, 
Eq. (51) has to be replaced in favor of 



£*,< — c 



Lr 1 



which is identical to the crossover scale in unconfined 
pulling. The longitudinal dynamics within OPT also 
changes. Note that neglecting the exponential in Eq. 
(88) is feasible only if the exponential approaches zero 
sufficiently fast, i.e. only if (ft) 1 / 2 ~ (fi*,<) 1/2 < 7 1/4 . 
For 7 < 7 C this condition is violated and we need to 



reevaluate Eq. ( 88 ) 



(Ap)optW 



d( l (, ,-2 9 2 ftA f I 2 







2 Vfi 

TT l p 



1 - e" 



q 4 + 7 q 4 + q 2 f + 7 



(93) 



where we neglected the second term in the integrand, 
which is small compared to the first one, given x <C 1. 
This is — again — identical to the unconfined case. In sum- 
mary confinement effects vanish once 7 falls below the 
critical confinement threshold j c . 



2. Release 

We will now show that the same conclusions carry over 
to the case of release. Due to large local values of the 
backbone tension, the presence of confining channel walls 
merely affects the dynamical characteristics of the release 
scenario within the OPT regime, where the projected 
end-to-end distance saturates at a constant value. More- 
over, just like in the case of pulling, sufficiently strong 
confinement shifts the crossover time scale t* below its 
unconfined counterpart such that the crossover to OPT 
(and — simultaneously — even complete relaxation of the 
chain) occurs at significantly earlier times. 

Again, all these conclusions are implicitly based on the 
assumption 7 3> j c - To see how confinement effects fade 
away for 7 <C j c> note that the validity of Eq. (60) is 
based on the assumption jt 3> 1, which is reasonable 
for t = 0(t t ) (£* given in Eq. (61)) provided 7 ^> j c . 



For 7 <C 7cj in contrast, we have to keep the neglected 
exponential in the second term [68 of Eq. ( 60 1 , which 
then becomes 




1 



-2q 4 t 



q 4 + i 

2 3/4 



(94) 



1/4 



Here we neglected 7 in the denominator of the sec- 
ond term, which affects the corresponding rational func- 
tion only for wave numbers q < 7 1 / 4 where the term 



(92) (l - e- 2 ? 4 *) 



J < ' ) is essentially zero. Moreover we neglected 

the first term in line two, which is of order 0(f~ 1 / 2 ) and 
therefore negligible compared to the second term, which 
is of order Oit 1 ! 4 ), for times t S> f~ 2 . Eq. (94) is again 
identical to the unconfined release scenario, whence also 
the time scale £*, which was determined in Eq. (61) on 
the basis of (Ap)opT, has to be replaced by 



C 4 L 8 



t*.< — 



(95) 



which is again identical to the one determined in Ref. 
[45] . In order to be consistent, this time scale has to be 
much smaller than 7 -1 , which again yields the condition 

7 < 7c- 



VI. CONTACT WITH EXPERIMENTS 

In order to make contact with experiments, this section 
presents the experimental time, force and length scales 
implied by our analytical findings. To this end we need 
to translate between effective confinement strength 7 and 
actual channel diameter D. Such a mapping can be per- 
formed by comparing the mean square transverse fluc- 
tuations of chains confined in a hard-walled channel to 
those trapped in harmonic confinements. This procedure 
yields 



7 



64a 4 p 



(96) 



where a « 0.17 has been determined by Yang et al. [2"5] 
by means of Monte Carlo simulations. Note that Eq. ( 96 ) 



is an equilibrium estimate, which we will extrapolate to 
non-equilibrium situations in what follows. 

We used this relation to exemplarily calculate typical 
time scales for DNA and F-actin molecules in prototypi- 
cal experimental setups for the three scenarios discussed 
in this work, cf. table |Vj Note, however, that the char- 
acteristic time scales strongly depend on the adjustable 
parameters L, 7 and f whence the respective time win- 
dows may be varied considerably in actual experiments. 
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Time Scale 


F-Actin 


DNA 


Pulling: 






r 2 


10~ 4 s 


10 _7 s 


V(7 1/2 f) 


0.01s 


10" 6 s 


A 


0.03 s 


0.09 s 


Release: 






r 2 


10- 4 s 


10~ 7 s 


A 


0.003 s 


0.02 s 


FRC: 






I" 1 


1.2s 


10~ 5 s 


A 


3.5s 


1.2s 



TABLE V: Characteristic time scales for the scenarios dis- 
cussed in this work. For F-actin chains we assumed £ p = 
16ttm, L = 13.5 ^m D = 2 fim and f = 2pN; for DNA we 
used £ p = 50 nm, L = 18.6 fxm D = 100 nm and f = 2 pN. 



Eq. 



(96) moreover allows us to translate the critical 
confinement strength 7 C into a critical channel diameter 
D c , above which confinement effects become invisible in 
the context of tension propagation and longitudinal re- 
sponse. Recalling j c = / 2 and using Eqs. (96) and (63) 
we find: 



D r = 2- 9 / 4 



c 



3/2 



£ 2 ' 



(97) 



For DNA molecules in typical experimental setups, D < 
D c is automatically fulfilled due to their small persis- 
tence lengths. Considering the much stiffer actin fila- 
ments (using L w 13.5 /jm and i p w 16 fim. |63j ) we find 
D c w 10 /xm, such that confinement effects are only visi- 
ble in rather narrow channels in this case. 

Finally, Eq. (96) may be utilized to investigate the 
accessibility of the limits x 1 and x <C 1. Note that 
/ 7 = K7 1 / 2 = k B T/2fi- 2 / 3 el /3 D- 4 ^, where / 7 gives 
the effective confinement force in original units. The 
forces applied in single molecule experiments are of or- 
der f ~ 1 pN [53] . Regarding F-actin as telling exam- 
ple (£ Pt actin ~ 16 /im [63] [65]), channel diameters be- 
tween D = 2 fim and D = 10 fim have been used in 
experiments on confined F-actin molecules |19j . giving 
rise to small force scale separation parameters of order 
x ~ 10~ 2 to i ~ 10~ 3 . On the other hand, channel di- 
ameters down to D = 30 nm [TS] and forces as small as 
f ~ 10 -2 pN [55] have been used in experiments on DNA. 
Using fj, DNA ~ 50 nm [T3EU we obtain x ~ 10 2 in this 
case. Both limits i« 1 and x ^> 1 are thus accessible in 
experiments. 



bending polymers. Our work was based on the investiga- 
tion of three paradigmatic scenarios including a sudden 
application and release of a point force f at the polymer's 
ends (confined pulling and confined release), as well as 
free relaxation from confinement (FRC). 

In the context of tension propagation and longitudinal 
response we showed that for both force scenarios dynam- 
ical signatures resulting from confinements become visi- 
ble only above a certain threshold for the effective con- 
finement strength 7, which is given by j c = £ 4 /(£ 4 L 8 ). 
This translates into an upper channel diameter (cf. Eq. 
(97)) which imposes almost no restrictions for experi- 
ments on DNA but requires channel widths below 10 /im 
for actin in order to observe effects caused by confine- 
ment. For weak confinement strengths the chain's dy- 
namics reduces to that of unconfined chains. For suffi- 
ciently strong confinement, we introduced the force scale 
separation parameter x — 7 1//2 /f in order to quantify the 
strength of confinement relative to the externally applied 
force f. We showed that both cases x ^> 1 and x <C 1 
are accessible in experiments and consequently derived 
the chain's dynamics for comparably strong [x <C 1) and 
weak (x ^> 1) forces f. The scaling laws governing the 
longitudinal dynamics of confined chains in the two force 
scenarios turned out to bear close resemblance with those 
of prestretched polymers [50 , suggesting to view confine- 
ment as effective prestretching force of strength ~ 7 1 / 2 
in this context. In particular, the polymer's longitudinal 
response (ARn)(t) saturates at its new equilibrium value 
already for times t > t* (cf. Eqs. (51) and (61)), even 
though t* is smaller than the usually obtained value for 
the longest relaxation time. This behavior is explained 
by a suppression of long wave lengths in the mode spec- 
trum of the stored length density due to confinement and 
has already been observed for prestretched chains where 
it was termed "premature saturation" |46j . 

The interpretation of confinement as effective pre- 
stretching force carries over to the scenario of free re- 
laxation from confinement (FRC) , which is for late times 
identical to the release scenario discussed in Ref. [45], 
safe for a different prefactor of (Ai?y) which arises as a 
consequence of somewhat different bulk conditions. For 
short times, however, the initially tension free contour 
builds up backbone tension in the bulk, which grows 
like ~ t 1 / 4 before it saturates at the constant value 
27 1 / 2 . This tension build-up in turn affects the longitudi- 
nal dynamics of the chain, which responds superlinearly 
«Ai?y) ~ i 9/s ) in this regime. The FRC scenario thus 
constitutes an intriguing example of initially straight con- 
tours, which were discussed quite generally in Ref. [STj] . 
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an ensemble average safe for subdominant contributions 



that depend on the boundary conditions for the contour 
[141 1151 ISO]. 

[68] The rational function in the first term varies on a scale 
of order 0(s/f), which is large compared to the scale on 
which the exponential approaches zero (0(t -1,/4 )), pro- 
vided t ^> f~ 2 , which we alway assume to hold (the op- 
posite case corresponds to small forces f < f c , where the 
chain reacts linearly to the externally applied force). 



